Brownian dynamics of a self-propelled particle in shear flow 
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Brownian dynamics of a self-propelled particle in linear shear flow is studied analytically by solving 
the Langevin equation and in simulation. The particle has a constant propagation speed along a 
fluctuating orientation and is additionally subjected to a constant torque. In two spatial dimensions, 
the mean trajectory and the mean square displacement (MSD) are calculated as functions of time t 
analytically. In general, the mean trajectories are cycloids that are modified by finite temperature 
effects. With regard to the MSD different regimes are identified where the MSD scales with t" with 
u — 0,1, 2, 3, 4. In particular, an accelerated (y = 4) motion emerges if the particle is self-propelled 
along the gradient direction of the shear flow. 

PACS numbers: 82.70.Dd, 05.40.Jc 



I. INTRODUCTION 

Mesoscopic colloidal particles which perform Brown- 
ian dynamics in a viscous solvent exhibit intriguing non- 
equilibrium behavior if they are exposed to shear flow. 
In fact, shear flow drastically affects their viscoelastic 
response [H-[l] , criticality and phase behavior [H-jl, and 
other collective effects with @, and without |ll Q3 
many-body hydrodynamic interactions. As an example, 
the mean square displacement of a Brownian particle 
in shear flow contains non-vanishing cross-correlations 
[l3l Il4j and shows several dynamical regimes when con- 
fined to an additional parabolic potential [TBj . 

Recently, active particles which are self-propelled by 
their own intrinsic motor (Ml have been studied. 
Even the single-particle dynamics is a non-equilibrium 
situation since the particles dissipate energy. Most of the 
recent studies of self-propelled particles were performed 
in a quiescent solvent (see, e.g., Refs. |19l - [25| h But also 
studies focused on self-propelled particles in an imposed 
shear field 26-29] gain more and more interest. Apart 
from its fundamental importance, the interplay between 
hydrodynamic flow and active particle dynamics is also 
relevant for ecology, having for instance the motion of 
bacteria in oceans, lakes, and rivers in mind [3fjj . Besides 
particles with a fixed shape, also swimmers which change 
their shape during propulsion were considered |3l( or a 
velocity field due to surface deformations was prescribed 
on the surface of the particle [HI, [H| . 

In this paper, we study the two-dimensional Brownian 
dynamics [34j of a single self-propelled spherical particle 
in linear shear flow based on the Langevin equations. The 
particle has a constant propagation speed which fluctu- 
ates in its direction and is exposed to a linear shear flow 
and a constant torque. Analytical results as functions 
of time t are presented for the mean trajectory and the 
mean square displacement (MSD). At zero temperature 
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the mean trajectories are cycloids, which are modified if 
finite temperature effects arc included. For the MSD dif- 
ferent regimes are identified where the MSD scales with 
t v with an exponent v = 0,1,2,3,4. While v = 3 is 
the most general case for long times, the other expo- 
nents are realized in the following special cases: v = 
in the absence of shear flow and fluctuations, v = 1 in 
the absence of shear flow, v = 2 in the absence of shear 
flow, torque, and fluctuations, v = 4 in the absence of 
fluctuations and for a special nonzero torque canceling 
the shear rotation. Live unicellular motile microalgae 
(Chlamydomonas Reinhardtii) that maintain their direc- 
tion in shear flow [35| provide an experimental realiza- 
tion of particles that resist the flow rotation for most 
of the time. The accelerated motion with the exponent 
v = 4 is found if a particle is self-propelled along the 
gradient direction of the shear flow. As a general result, 
the motion of self-propelled particles is greatly amplified 
by shear and, for special initial conditions, the constant 
acceleration can hugely enhance the mobility of individ- 
ual particles. Our results can in principle be verified in 
experiments of self-propelled colloidal particles in shear 
flow 0. 

The model considered in this paper has been solved 
before in two special limits. In the absence of shear, the 
solution was presented in Ref . [20j , while in the absence of 
self-propulsion, the traditional Brownian motion in shear 
flow is obtained [T3 . f37j | ■ It is important, however, that in 
none of these limits the exponent 4 emerges which is thus 
arising from a combination of self-propulsion and shear. 

In addition to the analytical analysis of mean trajec- 
tory and MSD, in this paper, we also focus on the com- 
plete probability distribution function for the displace- 
ment of a self-propelled particle in shear flow. As calcu- 
lated in simulation, it exhibits a transient double-peak 
structure, which is due to the self-propulsion and dis- 
torted by the shear flow. 

This paper is organized as follows: In Sec. |TT]we intro- 
duce our model of a self-propelled particle in shear flow. 
The analytical solutions are presented for the special case 
of zero temperature in Sec. Mil and for the general case of 



2 




©Me, 



FIG. 1. (Color online) Sketch of the self-propelled particle in 
Couette flow with shear rate 7. The system is considered in 
the two-dimensional (x, y) plane. A systematic torque M = 
A/e z and the direction u of the effective force F due to the 
self-propulsion are indicated. 

finite temperature in Sec. IIVI Section [V] contains some 
further results obtained by simulation and, finally, a con- 
clusion is given in Sec. I VII 



{g{t)9(t')) = 26(t-t')/(j3 2 D r ), where i,j e {x, y} and 
Sij is the Kronecker delta symbol. Here, we assume out 
of equilibrium noise to be close to the equilibrium one. 

The subsequent analysis in this paper shall be carried 
out in a dimensionless form. For that purpose, we in- 
troduce the following reduced parameters: x — x/R, 
y = y/R, t = D r t, = fc/F, g t = g/M, a = 4/3RF/3, 
fi = (3M, and £ = Pe r /2 = j/(2D r ), where Pe r = i/D r 
is the rotational Peclet number. With these parameters 
the new set of dimensionless Tangevin equations can be 
written as 

^ = 2£y + a(cos(0)+/ K ), (4) 
^ =a(sin(0)+;„), (5) 
g = - ? + M (l + p). (6) 

III. THE CASE OF ZERO TEMPERATURE 



II. THE MODEL 

In our model, we consider the completely overdampcd 
two-dimensional Brownian motion of a spherical self- 
propelled particle in Couette flow (see Fig. [T]). The dis- 
tance between the boundary plates and the particle is 
assumed to be large enough so that hydrodynamic in- 
teractions can be neglected. The mechanism of propul- 
sion is taken into account by means of an effective inter- 
nal force F = Fu — F (cos(<^>), sin(^))), which enters the 
Langevin equations. The centre-of-mass position vector 
r = (a;, y) and the angle <fi between the unit vector e x 
that is co-directional with the Cartesian x axis and the 
orientational unit vector u = (cos(0), sin(<^)) are used to 
describe the two-dimensional motion of the self-propelled 
particle. The effect of the shear flow enters the equations 
via the shear rate 7. Considering a systematic torque M 
as indicated in Fig. [I] a Gaussian white noise random 
force f = {fxify) and a Gaussian white noise random 
torque g leads to the following set of Langevin equations: 



Ax 

dt 
dy 
dt 



jy + PD t (Fooa(<f>) + f x 
PD t (FBm(4>)+f v ) . 



f t =-i+PD r (9 + M) 



(1) 
(2) 
(3) 



Here, f3 — l/(fceT) is the inverse effective thermal en- 
ergy and D t and D r are the translational and rota- 
tional short-time diffusion constants satisfying the ra- 
tio D t /D r = 4i? 2 /3 for spherical particles with radius 
R. With ( • ) denoting the noise average, the Gaussian 
white noise quantities are characterized by (fi(t)) = 0, 
(fi(t)W)) = 25ij8(t - t')/(/? 2 A), (g(t)) = 0, and 



At zero temperature, as there is no thermal motion, 
the noise terms in the system of dimensionless Langevin 
equations Q- © can be neglected and an analytical so- 
lution for the trajectory can be given. One has to distin- 
guish between the cases /i = £ and [i 7^ £. 

For \x ^ £ the reduced net torque u) given by to = ji — £ 
acts on the particle. This leads to the analytical results 



n , ( ~ acos(0n)" 
Ax = 2£ (y + l r 

+ ^ ~ ^ [sin(0 o ) - sin (<fo + ut)] , (7) 



Ay=— [cos(</>o) - cos (0 O + wr)] 

u> 



(8) 



where Ax — x(t) —xo and Ay — y(r) — yo with the initial 
reduced centre-of-mass position vector ?o = r(t = 0) = 
(xq, yo) and the initial orientation angle 4>o = 4>{t = 0) of 
the self-propelled particle. Equations an d © describe 
cycloids as can also be seen from the corresponding graph 
(dashed line) in Fig. f2J 
For the special case // = £ we get 

Ax = £a sin(0 o )T 2 + (2£j/o + a cos(0 o )) t , (9) 
Ay = asin(0 o )r , (10) 

Physically, the condition /j, = £ implies that the rota- 
tion of the self-propelled particle due to the shear flow 
is exactly compensated by the additional external or in- 
ternal torque represented by fi. Thus, the orientation of 
the particle remains constant all the time as there is no 
random torque in the case of zero temperature, either. 
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FIG. 2. (Color online) Mean trajectory of a self-propelled par- 
ticle in shear flow at finite temperature (solid line in the main 
figure) and the corresponding trajectory for T — (dashed 
line). The damping effect of the finite temperature is clearly 
visible. The parameters are a — 1 and <j>o — for all plots, 
while the normalized shear rate is £ = 5 in the main figure and 
£ = 25 in the upper inset. In the lower inset besides £ = 1, an 
additional torque determined by \i = 7 is considered as well. 



IV. GENERALIZATION TO THE CASE OF 
FINITE TEMPERATURE 

A. Mean trajectory 
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FIG. 3. (Color online) Mean square displacement of a self- 
propelled particle in linear shear flow. For all plots the 
strength of the self-propulsion is determined by a — 10 and 
the initial angle is <j>o = n/2. While T\ = for the dotted and 
the dashed-dotted lines, the solid and the dashed lines refer 
to situations with finite temperature T% > 0. A nonzero shear 
flow (£ = 5) is considered for the solid and the dashed-dotted 
lines. An additional torque (fi — 5) is taken into account for 
the dashed and the dashed-dotted lines. Depending on the 
temperature and the existence of shear flow and additional 
torque, different exponents of t can be identified in the long- 
time behavior. In the inset a close-up of the intermediate 
regime is shown. 



For finite temperature T > the Brownian motion due 
to collisions of the self-propelled particle with solvent par- 
ticles has to be taken into account. The mean trajectory 
of the self-propelled particle can then be written in the 
form 



(Ax) = 2£ 



(Ay) = 



l+LO 2 



1 + to 2 1 + LJ- 



1+UJ 2 



l+LO 2 



■e-\ (11) 



(12) 



Here, the parameters (f>Q and <fi^ containing the depen- 
dence on the initial angle 4>o were used. They are given 
by 



= sin(0 o ) + wcos(0 o ) , 
6(7 = cos(0 o ) - cjsin(0 o ) , 



<t>a =4> 



2£ 



1 



(13) 
(14) 

(15) 



The parameters A<p + , A0~, and Acj)^ are defined in 
exactly the same way, but with instead of 4>q in 
Eqs. (|13 [) -(|15 p . The analytical results for the mean tra- 
jectory are visualized by the solid lines in Fig. [5] and its 
insets. The main part of this figure illustrates basically 
two effects due to the finite temperature. Comparing the 



dashed line corresponding to T — with the solid line 
shows that the finite temperature leads first to a damp- 
ing in the amplitude and second to a reduced frequency 
of the oscillations. As presented in the lower inset, an 
additional torque leads to quite complicated mean tra- 
jectories. 



B. Mean square displacement 



Next, the MSD is discussed. The x and y components 
can be considered separately. The full analytical solu- 
tion, which is available, contains an enormous number 
of terms and is far too long to be given here explicitly. 
Therefore, we limit ourselves to presenting the integrals 
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and correlation functions that have to be calculated: 

<(A£) 2 ) = |eV+4^r 2 + ^ 

+ Aa 2 e fdn f dr 2 [dr 3 f dr 4 (sin {cp (r 2 )) sin (0 (r 4 ))} 
Jo Jo Jo Jo 

+ a 2 dn / dT 2 (cos (0 (n)) cos (<f> (r 2 ))) 
Jo Jo 

+ 8at?y T [ dn / dr 2 (sin (0 (t 2 ))) 
Jo Jo 



+ 4a£y T / dn (cos (n))) 
Jo 



+ 4a 2 £ / dn / dr 2 / dr 3 (sin (</> (r 2 )) cos (0 (r 3 ))) 
JO Jo Jo 



(16) 



Carrying out the calculations reveals a rich variety of 
regimes as visualized in Fig. [3] This is particularly re- 
markably because the underlying model situation is quite 
simple. In the following, situations that lead to the dif- 
ferent exponents v — 0,1, 2, 3, 4 shall be covered in more 
detail. 

The case v = is realized for T — in the ab- 
sence of shear flow. The corresponding motion of the 
self-propelled particle describes a closed circle due to a 
nonzero torque. The exponent v = 1 is the most gen- 
eral case without shear flow. Details of this model for a 
self-propelled particle were presented in Ref. 38]. The 
ballistic case v — 2 is more or less trivial and is realized 
for a zero torque at T = in the absence of shear flow. 
The particle simply moves on a straight line. The most 
comprehensive situation revealing the most interesting 
physics consists of systems with shear flow and torque at 
finite temperature T and in general leads to v = 3. In this 
case, basically three regimes can be identified. For short 
times, the MSD is linear in time as the simple diffusive 
motion is dominant. As soon as the self-propulsion be- 
comes significant, a crossover regime is observed, where 
contributions with t 2 , t 3 , and t are relevant. A more 
detailed discussion of this regime is given below. Finally, 
at the time scale 1/ 1 D r , the intermediate regime is termi- 
nated by the long-time t 3 law. While the exponent v = 3 
also occurs in the Taylor diffusion of passive particles [13] , 
the exponent v = 4, which is found for the situation with 
a torque exactly canceling the shear rotation also in the 
long-time behavior (see dashed-dotted line in Fig. 13]), is 
characteristic and only realized for active particles. 

In the inset of Fig. [3] a close-up of the crossover regime 
is shown. While the t 2 and t 3 terms in this transient 
regime also occur for self-propelled particles in a qui- 
escent solvent [381 ] . the t 4 contribution is only found for 
situations with shear flow. If the initial orientation of the 
particle is parallel to the y direction (</> = n/2) and the 
shear rotation is exactly compensated by the additional 
torque (/x = £), the motion is directed in the y direc- 
tion for a significant time. Thus, an accelerated motion 
with the exponent v = 4 arises. Unlike the case of zero 
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FIG. 4. (Color online) Time evolution of the probability dis- 
tribution function of a self-propelled particle. In plot (a) the 
distribution in the y direction is shown while in plot (b) the 
additional effect due to the shear flow in the x direction is 
illustrated. The different lines in blue, red, green, and black 
correspond tOT = l,T = 2,T = 3, and r = 6, respectively. 
The initial orientation of the particle points in the positive x 
direction. 



temperature, for nonzero noise the i 4 contribution is re- 
stricted to a transient regime, which is terminated by the 
general t 3 law. Although the condition /i = £ seems to 
be very special at first sight, an example of particles that 
resist the flow rotation for most of the time is given by 
live unicellular motile microalgae ( Chlamydomonas Rein- 
hardtii) that maintain their direction in shear flow [35| . 
Thus, the condition fi — £ is realized in experiment. 

The exact expression for the long-time behavior in the 
most general case (v — 3) is given by the first term in 
Eq. ([T6|) for passive particles. Depending on the strength 
of the self-propulsion and the effective torque in the shear 
flow, for self-propelled particles a second contribution be- 
comes relevant: 



hm ((Ai) 2 ) = ^ 2 



(l+W 2 ) 



(17) 



Going back to physical quantities, for the case of passive 
particles the leading term can be written as (2/3)7 2 £>ti 3 . 
Equation (|17|) verifies that the corresponding term for 
self-propelled particles is obtained by replacing the trans- 
lational diffusion constant D t for passive particles with 
the substantially enhanced long-time translational diffu- 
sion constant [19|, |38| for self-propelled particles with net 
torque oj. 



V. COMPUTER SIMULATION 

In addition to the analytical considerations presented 
so far, some further information can be gained by numer- 
ical computer simulation. While the different moments of 
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the probability distribution function for the displacement 
of a self-propelled particle in shear flow - in particular the 
mean position and the MSD - can be calculated analyti- 
cally, this is not the case for the probability distribution 
function itself. Thus, we make use of a Brownian dy- 
namics simulation here. Based on the Milstein method, 
a very good approximation for the probability distribu- 
tion function is obtained b y a veraging over a huge num- 
ber of particle trajectories [39| ■ Some results for the case 
that the rotational motion due to shear is exactly com- 
pensated by the additional torque are presented in Fig. S) 
The left plot (a) shows the distribution in the y direction, 
which is qualitatively different from the particle distribu- 
tion in the x direction [see Fig.0](b)] due to the fact that 
the consequences of the shear flow become obvious in the 
latter case. The double-peak structure, which is sym- 
metric in the left-hand plot, can be explained as follows: 
Given that the active particle is oriented in the x direc- 
tion at the beginning, no motion in the y direction due to 
the self-propulsion occurs at very short times. The im- 
portant aspect is that at finite temperature the particle 
does not only undergo translational, but also rotational 
Brownian motion. This results in a randomly changing 
particle orientation and thus gives rise to a varying y 
component of the propulsive motion. The rotational mo- 
tion of the particle at very short times due to collisions 
with the solvent molecules determines whether the self- 
propulsion drives the particle in the positive or negative 
y direction. As seen in the left plot in Fig. @J the re- 
sulting double-peak is clearly observed at times D r t < 3. 
For larger times, this detailed structure is more and more 
washed out and, finally, for very large times a Gaussian 
distribution function is established. Graphically speak- 
ing, this is the case when the initial change of the particle 
orientation becomes irrelevant due to the continuing ro- 
tational Brownian motion. As already discussed in Sec. 
IVI in this long-time regime the motion becomes diffu- 
sive again with a modified diffusion constant due to the 
self-propulsion of the particle. 

On top of the previous considerations, the effect of 
the shear flow becomes relevant with regard to the case 
shown in Fig. 2] (b). The form of the probability distri- 
bution function is not symmetric any more because the 
initial orientation of the particle points in the positive x 
direction and the shear flow induces an additional com- 
ponent of translational motion in positive x direction for 



positive y values and in negative x direction for negative 
y values of the current particle position. This explains 
the modified and asymmetrical form of the particle dis- 
tribution function in x direction. 



VI. CONCLUSION 

In conclusion, within an analytical solution of a two- 
dimensional model we have studied the mean trajectory 
and the scaling behavior of the mean square displace- 
ment (MSD) as a function of time if a self-propelled par- 
ticle is exposed to linear shear flow. For long times, the 
MSD scales with t 3 as a result of a combined diffusion 
and convection. If the self-progagation is initially ori- 
ented in the gradient direction of the shear flow, there is 
a constant-acceleration behavior where the MSD scales 
with i 4 , which finally crosses over to the i 3 law where the 
crossover time corresponds to the inverse rotational diffu- 
sion constant. The scaling behavior of the MSD and the 
transient double-peak structure of the probability dis- 
tribution function for the displacement of the particle 
as calculated in simulation can be tested in experiments 

For future work, it is interesting to generalize the 
present solution to the case of oscillatory shear 0, |4l|, l42j 
and to three spatial dimensions [38[ . Further future 
work should focus on the behavior of an ensemble of self- 
propelled particles which are coupled by both direct in- 
teractions and by hydrodynamic interactions mediated 
by the solvent [1, |43|, [44| . At finite densities, deviations 
from a Gaussian behavior [40| in the motion of parti- 
cles are expected due to the power-law decay of hydro- 
dynamic interactions [45j . The swimmer-tracer scatter- 
ing [45| of a self-propelled particle in shear flow might 
be interesting as well. An external imposed shear flow is 
expected to change the viscoelastic behavior drastically 
[29 . [27| and there is an intricate coupling between the 
intrinsic flow made by the motion of the active particles 
and the externally imposed shear flow. 
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